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Description 

Field of the invention 

s [0001] The present invention is in the field of digital radiography. The invention mor specifically relat s to a method 
of recognising one or more irradiation fields in a system wherein (a) radiation image(s) is (are) recorded on a photo- 
stimulable phosphor screen by means of an irradiation field stop (i.e. means for partially shielding the screen against 
irradiation). 

10 Description of prior art 

[0002] Radiologists may protect their subjects against unnecessary exposure to X-rays by use of X-ray opaque 
material. The material is placed in the path of the X-ray beam so as to shield those areas of the patient which are not 
regarded as diagnosticaily important. The resulting image contains relatively unexposed regions which comprise the 

is shadow cast by the X-ray opaque material. 

[0003] In addition, radiologists may use similar X-ray opaque material to protect whole sections of the resulting image 
from exposure. By subsequent transposition of the opaque material, they can achieve multiple, but spatially non-over- 
lapping, exposures of the same image. The motivation may be, for instance, to have two different views (hence sub- 
images) of the same subject juxtaposed on the same image. 

20 [0004] In computer-assisted radiology, the X-ray image is digitised as a two dimensional array of numbers, whose 
magnitude is related to the intensity of X-rays arriving at the detector, via the patient. The values in the array may be 
rescaled, prior to display or printing, in order to maximise the visual contrast in the area of diagnostic interest. Such 
rescaling depends, in part, on the intensity histogram of the X-ray image, or if applicable, its sub-images. To achieve 
this function automatically, best results are obtained if, for each sub-image, the shadow cast by the collimation material 

25 js excluded from consideration, since it cannot contain any useful information, but would otherwise dominate the lower 
end of the intensity histogram. 

[0005] The presence of the collimation shadow region can cause an additional problem in the display of the image 
on film. Therein, the shadow region is relatively bright, and if displayed unmodified, may impair diagnosis of subtle 
lesions due to dazzle, especially if the unexposed region is relatively large. 
30 [0006] The following earlier work refers only to the situation where the whole X-ray image is subject to a single 
exposure (i.e. not multiply exposed). 

[0007] In U.S. patent 4,g52,807 by Adacht a method is disclosed for adjusting the image processing conditions by 
selecting only those pixels which correspond to the object or the background, i.e. the 'signal 1 region, under the assump- 
tion that the collimation material does not partly cover the object, i.e. the shadow region does not touch the object of 
35 the image, but rather its surrounding background. In many kinds of examinations however, parts of the patient will be 
shielded from exposure to avoid useless X-ray load, as is the case with lumbar spine and shoulder examinations. 
[0008] U.S. patent 4,804,842 discloses a method for removing the lower part of the histogram corresponding to the 
pixels in the shadow region in order to adjust the image processing conditions. 

[0009] However, techniques for histogram modification atone, do not address the second problem, i.e. that of dazzle. 
40 its solution demands an explicit detection of the shadow region, as distinct from the signal region. The offending shadow 
region may then be excluded from the display process. In addition, the first problem of rescaling for optimum visual 
contrast can then proceed on the signal region alone, without the complicating influence of the shadow region on the 
intensity histogram. 

[0010] Previous work in the field of automatically delineating the signal/shadow regions has relied on piecemeal 

45 attempts at the solution. 

[001 1] Many of the proposed methods are based on collecting candidate signal/shadow border pixels (i.e. edge pixels 
at the interface between signal and shadow regions) by thresholding on first differences along predetermined lines that 
might cross the signal/shadow border. A typical example of this approach is proposed by Funahashi in the U.S. patent 
4,970,393. The signal region is recognised as a region enclosed by pixels where the first difference exceeds some 

50 specified threshold, either positive or negative. 

[0012] A refinement of this technique is proposed in European patent 0 342 379, where multiple candidate pixels 
along a searching line are ranked, and only the highest rank candidates on each searching line are consisered in 
assembling the signal/shadow boundary. In both methodes the underlying assumption is that the signal/shadow bound- 
ary has strong contrast everywhere. In many cases however, contrast along parts of the signal/shadow boundary is 

55 nearly zero, especially if the collimation material is partially covering very dense parts of the body. So if the recognised 
signal/shadow boundary is assembled from high difference pixels many gaps will occur in general. 
[0013] In European patent 0 285 174 a straightforward application of the basic Hough transform technique is pro- 
posed to delineate the signal/shadow boundary, which is assumed to be polygonal. Th Hough transform is applied to 
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prospective edge points obtained by differentiation. The unterlying assumptions ar that th signal/shadow boundary 
is piecewise linear, and has high contrast with respect to the contrast of line features within the signal r gion. Howev r, 
experimental evidenc shows that edges within th signal region, e.g. bone edges, frequ ntly have high r contrast 
than the signal/shadow boundary. As a consequence points in Hough transform space that correspond to straight bone 
5 boundaries may have higher scores than some of th actual signal/shadow borders, in which case the subsequent 
delineation procedure will fail. 

[0014] A different approach is disclosed in the European patent 0 360 231. The image is subdivided into mutually 
adjacent non-overlapping tiles. A statistical parameter is computed indicating the edge content of the tile, either by 
differentiation in image space, or in frequency space. The tiles are accordingly classified as belonging to the relevant 
10 image region or not. 

[0015] We ran elaborate experiments and concluded that single criteria based on edge contrast, or local statistical 
parameters at both sides of a prospective signal/shadow boundary, such as mean value, minimum, maximum, variance 
or spatial frequency coefficients are inadequate in delineating the irradiated field with a high level of confidence, say 
> 98 percent. We found that the signal/shadow border contrast will vanish in some cases; anatomical structures may 

is be straight and have high contrast edges; the image signal (which represents local dose) in the shadow region will 
occasionally be higher than the signal in the unshielded (i.e. diagnostic) image region - in cases where the scattered 
radiation is high and some parts of the body are difficult to penetrate by radiation. Also we found that there is no 
essential distinction between image areas corresponding to very poorly penetrated bone, and areas conesponding to 
shielded parts. Both kinds of areas will be characterised by an extremely bw mean value (i.e. dose), similar noise level 

^0 (both quantum limited), and by the absence of other features (since the signal to noise ratio is very tow). 

[001 6] I n order to obtain a high success rate, a preferred method should consider many hypotheses as to the location 
of the signal/shadow boundary (also called prospective signal/shadow boundaries), and keep or reject each hypothesis 
"as a whole", instead of keeping or rejecting individual components of a hypothesis, such as edges. 
[0017] Multiple hypotheses as to the location of the signal/shadow boundary are considered in the method disclosed 

25 jn U.S. patent 4,952,805. For each of the corresponding candidate signal regions, two signal histograms are computed : 
one for the pixels inside the candidate irradiation field and one for the pixels outside. The candidate boundary which 
is characterised by the largest degree of interclasses separation as determined from both histograms is selected to be 
the recognised irradiation field. The degree of interclass separation is defined in such a way that a high value will result 
if there is very little overlap between both histograms. Hence this criterion will be effective only if the range of signal 

30 values within irradiation field are well separated from the signal values of the shielded region, which does not hold in 
many practical cases as we verified on the basis of digital radiographs from various examinations. Also, there is a 
practical limitation regarding the feasibility of the above mentioned approach, in that the number of hypotheses as to 
the location of the irradiation field may be prohibitively large. 

[0018] Previous work in general has focused on methods in which finding the boundary of an irradiation field essen- 
35 tially relies on a very small number of mostly local criteria Reviewing a wide variety of examinations confirms that 
individual criteria will only work in a very limited number of cases. A high success rate can only be achieved if many 
sources of evidence are combined, ranging from local image features to global constraints. 

[0019] According to European patent application EP-A-610.605, filed by Agfa-Gevaert N.V., many hypotheses as to 
the location of the boundary are generated and the final decision as to the correct hypothesis is contingent on the 
40 results of a number of high-level tests applied to each. 

[0020] EP-A-610.605 discloses a method for the determination of the location of the signal/shadow boundary in an 
X-ray image represented in a digital signal representation that comprises the steps of : 

i) Extracting low-level primitives from the X-ray image {X(i,j)}» 
*s jj) Forming a reduced number of intermediate-level primitives from the low-level primitives, 

iii) Building hypotheses as to the location of the signal-shadow boundary from combinations of intermediate-level 
primitives, during which each combination is subject to intermediate-level tests, whereupon partial or complete 
hypotheses are rejected or accepted, 

iv) Performing high-level verification tests on each hypothesis, whereupon hypotheses are rejected, or accepted 
50 at some cost, 

v) Selecting the hypothesis with the least cost. 

[0021] The behaviour of the radiologist, in placing the collimation material around the patient, and in contriving to 
procure multiple exposures, may broadly be characterised by a set of simple rules. This invention addresses only the 
55 use of collimation material which casts a shadow whose outline is piece-wise linear. More specifically, for each individual 
exposure, the boundary between signal and shadow regions is assumed to be a rectangle of arbitrary location and 
orientation. It is assumed that the signal region is on the interior, and that any part of the rectangular boundary may, 
or may not, have a protection onto the (rectangular) detector, and therefore the image array; see figure 1. In addition, 
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it is assumed that the piac ment of X-ray opaque material (masks) to permit multiple exposure is always to divide the 
image into two, approximately equal, parts, either vertically or horizontally aligned; see figure 2. Not that the combi- 
nation of multipl -exposur , and patient shielding (collimation) can giv a wide variety of xposure area shapes; see 
figure 3. 

5 

Objects of the invention 

[0022] It is an object of the present invention to provide a method of recognising multiple, independently acquired, 
irradiation fields in a single radiation image, in a system wherein a radiation image has been multiply exposed using 

10 irradiation field stops as masks. 

[0023] It is a further object to provide a method for determination of the location of the multiple-exposure boundary 
of an X-ray image, wherein the X-ray image is represented by a digital signal representation and wherein the multiple- 
exposure boundary is the locus of the interface between the two or more sub-images resulting from the use of masks. 
[0024] It is a still further object to provide such a method for application in a digital radiography system. 

is [0025] Still further objects will become apparent from the description hereafter. 

Statement of the invention 

[0026] To achieve the above objects the present invention provides a method of determining a multi-exposure bound- 
20 ary between sub-images in a multi-exposure X-ray image represented by a digital signal representation and of deter- 
mining a signal/shadow boundary in each sub-image, comprising the steps of 

i) Extracting low-level line primitives from the X-ray image, 

ii) Forming intermediate-level primitives from the low-level primitives, 

25 iii) Building hypotheses as to the location of the multiple-exposure boundary from combinations of intermediate- 

level primitives, 

iv) Associating with each hypothesis a cost, 

v) Subjecting each hypothesis to evaluation tests via a rule-based partitioning reasoning network, 

vi) Selecting the hypothesis with the least cost, 

30 vii) Partitioning the image into sub-images according to the selected hypothesis, and 

viii) Submitting each sub-image for subsequent detection of the signal/shadow boundary to the following steps (ix) 
to (xv): 

ix) Extracting low-level primitives from each sub-image, 

x) Forming intermediate-level primitives from the extracted tow-level primitives, 

35 xi) Building hypotheses as to the location of the signal/shadow boundary from combinations of intermediate-level 

primitives, 

xii) associating a cost with each hypothesis, 

xiii) Subjecting each hypothesis to evaluation tests via a rule-based collimation reasoning network, 

xiv) Selecting a hypothesis with the least cost, 

40 xv) Construction of a binary-valued mask image representing signal and shadow area to be used to extract the 

signal area in the X-ray image. 

[0027] The invention includes two reasoning networks (rule-based systems) which, at each node, contain the rules 
for construction, costing, and evaluation of hypotheses of both multiple-exposure and signal-shadow boundaries. A 
45 different set of rules is used for detecting the partitioning boundary and for detecting the signal-shadow boundary in a 
sub-image. The first set of rules is called the partitioning rule network, the second is the collimation rule network. 
These networks have the following properties: 

a) Arrangement of the reasoning network so as to first classify hypotheses as topological, and subsequent as 
50 geometrical, archetypes. 

b) Arrangement of network to subsequently further refine and evaluate hypotheses so classified. 

c) Ability to train the network, by changing rule behaviour interactively, and monitor consequences of that change 
on the subsequent evaluation of the current hypothesis. 

55 [0028] The main departures from the method disclosed in European patent application 610 605 are : 

A) The simultaneous detection of multiple-exposure and (possibly more than one) signal/shadow boundaries in 
each X-ray image. 
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B) The use of reasoning networks embodying archetype classification, refinement, and evaluation of the hypoth- 
eses for the location of multiple-exposure and signal/shadow boundaries. 

[0029] Th method permits codification of xpert-syst m type knowledge as to the geom try, topology, and statistics 
5 of the multiple-exposure and signal/shadow boundaries, and associated regions. 

[0030] The invention is for example, applicable to a digital image signal obtained by exposing a photostimulabfe 
phosphor screen to an image originating from a partially shielded exposure of an object to X-rays, by scanning the 
screen with stimulating irradiation, detecting the light emitted upon stimulation and converting the detected light into a 
digital signal representation. 

10 [0031] Commonly the digital signal is subjected to image processing such as contrast enhancement and signal-to- 
density mapping before being applied to an output device for reproduction or display. The result of the process of 
contrast-enhancement of an X-ray image is improved by permitting such to perform separate calculations on each of 
the diagnostics ly useful signal parts of each of the sub-images of the X-ray image. The results of the separate calcu- 
lations may then be applied to each sub-image independently to present the signal part of the image with an optimal 

15 contrast, and also to reduce the brightness in shadow region so as to avoid dazzle during diagnosis. 

[0032] In the statement of the invention and in the description hereinafter the term "X-ray image' is used in the 
meaning of a two-dimensional digital image representation as an array of numbers, the magnitudes of which are related 
to the intensity of the X-rays arriving at a 2-dimensional physical detector such as a photostimulable phosphor screen. 
[0033] By "primitives* is meant features, such as lines (see figure 5), and corners, intersection points etc. (see figure 

20 6) which are extracted from the image as the result of low-level image processing. 

[0034] An "intermediate-level primitive" is a combination of (basic) primitives which may be considered collectively 
when building a hypothesis. Intermediate-level primitives are generally fewer, and (spatially) larger than are (basic) 
primitives. Each primitive is constructed from a mutually exclusive set of pixels. Each intermediate-level primitive is 
constructed from a mutually exclusive set of (basic) primitives. As an example, a combination of lines (primitives) into 

25 line-clusters (intermediate level primitives) is illustrated in figure 7. Here, the line-cluster (figure 7 (c)) is the result 
of two lines l 3 and l 4 (figure 7 (b)), which have been grouped together into a single entity. 

[0035] By the term "signal region" is meant the spatial area of an image which is not obscured by collimation material 
from a direct line-of-sight of the X-ray source, see figure 1b. 

[0036] The "shadow region" is the area of the image which is obscured by the collimation material from a direct line- 

30 of-sight of the X-ray source; see also figure 1 b. 

[0037] The "multiple-exposure boundary* is the projection of the edge of the X-ray masking material on the image 
(see figure 2). There may be cases where there is little or no contrast across that projection, as a result of collimation 
material in each of the sub-images, as in figure 3 (a) and figure 3 (c). In that case, by "multiple-exposure boundary" is 
meant any straight line which splits the image into two sub-images, without crossing the signal region in either sub- 

35 image. 

[0038] The "signal-shadow boundary" is the locus of the interface between the signal and the shadow regions. E.g. 
the pixels at the interface between the signal and the shadow regions. 

[0039] Steps i) and ii) of the method of the invention achieve a reduced set of primitives which can be assembled 
into hypotheses as to the location of the multiple-exposure boundary. Similarly, for each sub-image, steps vi) and vii) 
40 of the method of the invention achieve a reduced set of primitives which can be assembled into hypotheses as to the 
location of the signal/shadow boundary. 

[0040] In the context of collimation detection, a "hypothesis" is defined as a possible solution proposed by the method 
of the present invention for the segmentation (division into similar regions) of the original image into signal and shadow 
regions, e.g. a two-dimensional array of numbers isomorphic with the X-ray image, with pixel labels for signal and 
45 shadow regions. 

[0041] Hypotheses are generated on the reduced (intermediate-level) set because, in general, the number of low- 
level primitives is too numerous to perform an exhaustive evaluation (steps iv and x) of every possible combination 
(steps iii and ix). 

[0042] With reference to steps v) and xi) : following an exhaustive search (including evaluation) by the reasoning 
50 network of all topological combinations of intermediate-level primitives, the hypothesis with the least cost is the one 
adopted as the solution. 

[0043] At step xii), by combining the solutions for the multiple-exposure boundary, (if applicable), and solutions for 
the signal/shadow boundaries (in each of the sub-images), a final result is constructed which is a mask image with 
(integer) labels at each pixel which is unique to (each of) the signal region(s), and a label at each pixel which is unique 
55 to the shadow region(s). 

[0044] By "reasoning network" at point a) is meant a directed graph wh rein each node corresponds to a binary- 
valued rule which acts upon the current hypothesis, and wherein each directed arc, from a parent to a child node say, 
is the flow of control for the current hypothesis from parent to child, provided the rule-value at the parent is true - see 
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figur 8. 

[0045] There are two, effectively independent reasoning networks, each with its own entry point and terminal nodes. 
The two networks are lor the evaluation of multiple-exposure hypothes s, and for the evaluation of collimation (signal/ 
shadow) hypotheses. 

s [0046] By "archetype" at point a) is meant (a hypothesis belonging to) one of a set of generic - topological and 
geometrical - configurations, which may be subjected to an associated set of rules corresponding to a major branch 
in the reasoning network. As implied at point b) : within each of these major (archetypical) branches, the rules are 
designed to confirm, refine, and evaluate the archetypical classification, insisting on particular geometrical and statis- 
tical properties of the hypothesis known to be associated with that archetype - see figure 9. 

10 

Brief description of the drawings 

[0047] Particular aspects of the present invention as well as preferred embodiments thereof will be illustrated by 
means of the following drawings in which 

75 

Figure 1 shows signal and shadow regions for typical single-exposure image, 

Figure 2 shows examples of multiple exposures, 

Figure 3 shows examples of multiple exposure and collimation combined, 

Figure 4 shows the placement of multiple-exposure boundaries, 
20 Figure 5 illustrates the definition of edges, 

Figure 6 shows low-level primitives, 

Figure 7 illustrates line-clustering. 

Figure 8 is an example of a reasoning network, 

Figure 9 are instances of signal/shadow boundary archetypes, 
25 Figure 10 illustrates the derivation of multiple-exposure boundary configurations, 

Figure 11 illustrates the relation between line-clusters, multiple-exposure boundary search tree, and hypotheses, 

Figure 12 illustrates the relation between line-clusters, signal/shadow boundary search tree, and hypotheses, 

Figure 13 shows a complete set of archetypes of the multiple-exposure reasoning network, 

Figure 14 shows a complete set of archetypes of the collimation (signal/shadow) reasoning network. 

30 

Detailed description: 

Steps i) and ii) Line detection and clustering 

35 [0048] These steps have already been described in detail (also as steps i and ii) in European patent application EP- 
A-610 605. 

[0049] The steps are unchanged from the system described in EP-A-61 0 605. They form the low-level image process- 
ing and analysis applied to the whole image. 

[0050] In this system, however, the results are then applied to the problem of detection and localisation of the multiple- 
40 exposure boundary (steps iii, iv and v), whereas in the earYier system they were applied to the detection and localisation 
of the signal/shadow boundary (here : steps ix, x and xi). 
[0051] Steps (i) and (ii) are as follows: 

Step i) Extraction of low-level primitives 

[0052] The preferred low-level primitives are lines. They are derived from the X-ray image in two stages: edge- 
detection and segmentation. 

1.1 Edge Detection 

50 

[0053] The purpose of this stage is to produce an edge-image. The system according to the present invention incor- 
porates a Canny edge-detection process. 

[0054] The following are definitions of some terms used in connection with a Canny edge detection process: 
[0055] Canny edge-detection: A particular process for edge-detection comprising optional smoothing, followed by 
55 production of an edge-image wherein a pixel is labelled as an edge-pixel if and only if the Sobel squared magnitude 
at that point is a local maximum along the (one-dimensional) trajectory which is normal to the Sobel orientation. 
[0056] A Canny edge-image (E(r)} is formed by non-maximal suppression of the Sobel squared magnitude in the 
direction of th orientation of the image gradient, i. . normal to the Sobel orientation: 
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E(r) = 1lfG(r)>G(r + R q )^G(r-R q ) 
= 0 otherwis 

5 

where 

R q = (1,0) if IIQ(r)l - n/2 ISrc/8 
W = (1 ,-1 ) if IIQ(r)+rc/4l - n/2 I < n/8 

= (1,1) if IQ(r)+K/4kn/8 
1S =(0,1) if IQ(r)l<;*/8 

[0057] Let the vectors r, r q represent the two-dimensional coordinates on the image plane (r=(i j) say). 
[0058] Smoothed image: An intrinsic image derived from the X-ray image wherein the high spatial frequencies are 
suppressed relative to the low spatial frequencies. E.g. the result of convolution of the X-ray image with a (two-dimen- 
20 sional) Gaussian. 

[0059] Sobel squared magnitude: The pixel-value G(i,j) representing the square of the gradient-magnitude of the 
image {X(i,j)} as determined by the following formula: 

2S = (X(i+1,j) - X(i-t,j)) 2 + (X(i,j + 1) - X(i,j-1)) 2 

[0060] Sobel orientation: The pixel-value Q(i,j) representing the orientation of the tangent to the iso-gradient-mag- 
nitude of the image {X(i,j)}, as determined by the folbwing formula: 

Q(i,j) = atan2(X(i+1 ,j) - X(i-1 ,j), X(i,j+1 ) - X(i,j-1 )) 

where atan2(y,x) is the arctangent C language function defined in [-jc,n]. 

[0061] The edge-detection process incorporated in our system comprises four stages: 

35 

1.1.1 Subsampling of the X-ray image 

[0062] The subsampling reduces the number of pixels of an image. For example: an original image obtained by 
reading an irradiation image stored in a photostimulable phosphor screen comprising approximately 5 million pixels, 
40 is reduced to an image of about 20000 1 2-bit pixels, preserving the aspect ratio of the original image. 

1 .1 .2 Canny Edge-Detection 

[0063] The subsampled image is passed to a Canny edge-detector to produce images: Sobel-squared magnitude 
45 Sobel-orientation {Q(i,j)}, and edge-image {E(i,j)}. The Gaussian smoothing is via a symmetric 5x5 convolution 

mask with weights calculated to effect a standard deviation of 0.35 pixels. 

1.1.3 Elimination of Weak Edges 

so [0064] Edge-points with a low gradient magnitude indicate a weak edge. Studies have shown that the signal/shadow 
boundary is generally of high contrast, giving a large gradient magnitude. Very weak edges are therefore unlikely to 
form part of the signal/shadow boundary. Those with a Sobel-squared-magnitude below a threshold of 2000 are re- 
moved from {E(i,j)} f and correspondingly from {Q(i j)}, and {G(iJ)}. 

55 1.1.4 Angle of Locus 

[0065] The orientation {0(i,j)}, of the locus of edge-points in the edge-image {E(i,j)}, is calculated in three stages: 
edge-thinning, connectivity analysis, linear regression. 
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1.1.4.1 Thinning 

[0066] The edge-image is destruct'rv ly thinned. The motivation is thr efold: 

5 1 . The parameters of a regressed straight line, including the goodn ss-of-fit, are more accurately determined from 

the spatial distribution of edge-points. 

2. The various types of junction (TERMINAL, T-JUNCTION, ...) are more easily detected. This is required for corner- 
detection (below). 

3. There then exists a close correspondence between the length of a straight line, and the number of (connected) 
10 edge-pixels of which ft is made up. This is required for line-statistics (above). 

[0067] Thinning is achieved by removing edge-pixels which do not affect the local (3x3) connectivity of other edge- 
pixels. 

is 1.1 .4.2 Connectivity Analysis 

[0068] Connected edges are uniquely labelled by recursive edge-walking to create a label image {L(i,j)}. All edge- 
pixels mutually connected share a common, but otherwise unique, label. 

20 1.1 .4.3 Linear regression 

[0069] In this connection the following definitions are first given: 

Linear regression: The process by which an edge can be optimally parameterised according to the spatial location 
2S of its constituent edge-pixels. E.g. the determination of the parameters of a line which are the radius to the normal 

(from the origin), and its orientation relative to the X-axis, by a least squares (of perpendicular distances) fit. 

Weighted linear regression: The process by which an edge can be optimally parameterised according to both the 
spatial location of its constituent edge-pixels, and the degree of certainty as to their correct classification (labelling). 
so E.g. the determination of the parameters of a line which are the radius to the normal (from the origin), and 

its orientation O^r relative to the x-axis, by the minimisation of a weighted least squares (of perpendicular dis- 
tances) fit, the weighting being the Sobel squared magnitude. 
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[0070] Let p k be the positions of the edge-pixels (indexed by k) which comprise the line: 



n , rt M n <yh 
Pk^Pk 'Pk > 



and let the positional covariance matrix be n 



r _ p(x,x) r (x,y) 

r <y*) r (y.y) 

given by I^>^k (a) Pk (b) Wk " ZkPk (a X^k^ )w kfor In then Qw LR = 1/2 atan2(-2r< x 'y>, r^-iiy.y)) and F^ LFl 
= (sinfQyyLRjS^Mwk + c^Qw^Xtfyto)^)/!^* The eigenvalues \+K of T are the sum of square errors (i.e. pro- 
portional to the variances) in each of the principal directions ie along, and normal to, the best fit straight line: 

2X. =Tr(T)± (Tr(r) 2 -4Det(r)) 1/2 



[0071] In the described embodiment of the present invention the orientation of the locus at each point is calculated 
by linear regression of connected edge-points in a local neighbourhood (nominally 7x7). The connectivity is guaranteed 
ss by including only those pixels which, in {L(i,j)}, share a common label. The positional-covariance matrix gives the 
orientation of the maximum varlanc of the data which is also the orientation of the least-squares fit straight line with 
two degrees of freedom. 

[0072] With reference to the definition of weighted linear regression, with w^l : 0(iJ)=Q WLR is th orientation of the 
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locus of edge-points without regard to the sense of the original edge. The ambiguity in the sense of the locus is resolved 
by borrowing from th Sobel orientation Q(r): the direction of the tangent to th Sobel squared magnitude iso-contour 
is always close to the direction of th tangent to the locus of edge-points as determined by the Canny edge-detector. 
Therefore the orientation is adjusted to minimise the angular separation between 0(i,j) and Q(i,j): 

5 

Min (O(i.j) - Q(i.j) + kn> -> kmin 

k in {-3. .3) 

10 0<i.j) :- O(i.j) - Q(i.j) + KLin* + 2h * 

where n is chosen so that the result is in [0,2k]. 
1 .2 Segmentation 

15 

[0073] Line primitives are built from the remaining edges. The chains of edge-pixels are broken into piece-wise 
straight sections, and the pixels from each section assembled into a line structure. There are four steps to this proce- 
dure: 

20 1.2.1 Vertex Analysis 

[0074] An intrinsic image {V(i,j)} is constructed with labels for the classification of edge-pixels in {E(ij)}. According 
to the 3x3 connectivity of edge-pixels in (E(i,j)}, one of the following labels are awarded to the corresponding point in 
{V(i,j)}: ISOLATED-POINT, TERMINAL, EDGE-POINT, T-JUNCTION, X-JUNCTION. 

25 

1 .2.2 Comer Detection 

[0075] Corners are detected by a refinement of the EDGE-POINT labels in {V(i,j)}. There are two stages: 

30 1 .2.2.1 Curve detection 

[0076] An intrinsic image {U(i,j)} is calculated with values assigned to each pixel which corresponds to an EDGE - 
POINT labelling in {V(i,j)}. The pixel values in {U(i,j)} measure the local conformity of edge-pixels in {E(i.j)}, in a 7x7 
neighbourhood, to a straight-line. The positional covariance matrix is calculated as described in 1 .1 .4.3 with all w^=1 . 
35 [0077] The conformity, (to a straight line) is taken to be the ratio of major to minor eigenvalues of the positional 
covariance matrix, ie \ + fk_. Pixels previously labelled as EDGE-POINTs in {V(i,j)} are then relabelled as CURVEs, if 
the conformity is less than a threshold value: 
If V(i,j)=EDGE-POINT , and \X+ > 0.04, then V(i,j) = CURVE. 

40 1.2.2.2 Non -maximal supression 

[0078] The pixels in {U(i,j)} are identified as comers if they are both CURVE points in {V(i,j)} (i.e. having a value 
above some threshold), and they are locally maximum in a 3x3 neighbourhood. The corresponding pixels in {V(i,j)} are 
relabelled as CORNERS. 

45 

1 .2.2 Extraction of connected edges 

[0079] Those pixels in (V(i,j)} still retaining an EDGE-POINT classification are candidates for incorporation into the 
low-level, primitives. Such pixels are extracted from {V(i,j)}, into (S(i,j)} say, which is the linear segmentation as an 
50 intrinsic image: each connected chain of pixels is given a unique label in {S(i,j)J. Each connected set of pixels is then 
gathered into a line structure which is the low-level primitive. 

1.2.2.1 Rejection of short lines 

55 [0080] Studies have shown that sides of the signal/shadow boundary generally have at least some portion of their 
length for which the edge-detector furnishes an unbroken chain of edge-pixels. In order to avoid false positives there- 
fore, an implicit rule is that a side cannot be composed entirely of very short edge-chains, but must have at least one 
line (in the line-group) which is of substantial length. Prior to segmentation and extraction of the line primitives therefore, 
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chains of edge-pixels shorter than some threshold length are rejected (from {S(i j)}) at this point. By this means, th 
number of low-level primitiv s is reduced without compromising the accuracy: short lines may form part of th final 
hypothesis foil wing data-retrieval. 

s 1 .2.3 Weighted linear regression 

[0081] A weighted linear regression is performed on each of the lines, by which means the optimal parameters 
(normal radius to origin, and orientation), and conformity are recorded. The weights are the Sobel-squared magnitudes 
of the contributing edge-pixels, as a fraction of the maximum Sobel-squared magnitude: 

10 

v k (i.j) - G(i.j)''&ax lU(u.v)} V i.j 

v u.v 

is 

1 .2.3.1 Rejection of poor lines 

[0082] Lines are no longer 'local' entities. The test on conformity to a straight line.is now repeated with a less tolerant 
threshold in accordance with the extra precision afforded by the (arbitrarily large) number of edge-pixels in the line: 

20 

1 .2.4 Linear extrapolation 

[0083] The surviving lines are extrapolated from each end, until they intersect either the edge of the image, or an 
edge-point which does not belong to that line (as determined from the labels in {S(i,j)}). The 'limits' so determined, are 
2S saved in the line structure. A distinction is made between each of the two limits which correspond to the two ends of 
the line. This distinction respects the sense of the line. 

1.2.5 Data-retrieval 

30 [0084] Edge-points are now recovered from {E(i,j)} which were not classified as EDGE-POINTs in (V(i,j)} (perhaps 
because they were reclassified as CORNERS, CURVES or JUNCTIONS etc). For each line, two rectangular areas are 
scanned in (E(i,j)} each of which is centred on (aligned with) the projection of the line as given by the regressed pa- 
rameters RyY LR and O^r. The rectangles extend from the line end-points, to the limits as determined by extrapolation. 
The width of the rectangle is 5 pixels. All edge-points in {E(i,j)J which were not already attributed to the line are now 

35 included. 

Step ii) Extraction of intermediate-level primitives 
2.1 Linear agglomerative clustering 

40 

[0085] In this connection the following definition are given: 

[0086] "Linear agglomerative clustering 1 generally is an iterative process for grouping together lines which are 
deemed (by some criteria) to be sufficiently similar. E.g. a pair-wise agglomerative clustering process based on an 
analysis of variance test, and broadly as described in "Unsupervised Estimation of Polynomial Approximations To 
45 Smooth Surfaces in images or Range Data" by J.F. Silverman and D.B. Cooper, Proc. IEEE Cont. on Automation and 
Robotics, 1986, page 299-304. and defined hereinbelow, but with the mergable entities being lines rather than regions, 
and the constraint that two mergable regions be touching replaced by a constraint that two lines (or line-clusters) be 
linearly coherent. 

[0087] "Analysis of variance" is a statistical method for comparing (the similarity of) two sample distributions under 
50 the.assumption that the underlying distribution is Gaussian; but the variance is unknown. The method determines an 
F-statistic, which is the ratio of the two sample variances, correctly normalised to take account of the true degrees of 
freedom. 

[0088] When the test is to determine the likelihood of two distributions having the same parent distribution, then the 
ratio is between the unexplained and the explained variances. When applied to pair-wise clustering, the unexplained 
55 variance is the extra that is incurred by the proposed merge: 

<°u 2 > = (X m 2 -X a 2 'X b 2 V|v m -v a -vJ 
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where a,b and m denote the separate, and merged entities respectively; x and v are the chi-squared variates and 
corresponding degrees of freedom. The explained variance is that prior to the proposed merg : 

5 «^e 2 > = (Xa 2 + X b 2 )/(v a + V b ) 



10 



15 



20 



and 



F ( v o- v e) = <a u 2>/< °e 2> 



is the F statistic, having degrees of freedom 

v u = | v m- v a- v d 

and 



[0089] The F statistic is used to index the cumulative distribution function Q(Frv a ,v b ) which gives the probability of 
that value or greater occuring by chance. If the unexplained variance (as a fraction of the explained) is low then F is 
small, and Q(Flv a ,v b J) is large. In that case, the entities may be safely clustered. 
25 [0090] In the described embodiment of the present invention lines undergo linear agglomerative clustering according 
to a test of their linear coherency. For all line or line-cluster pairs, if the pair passes the test of linear coherency, an F 
statistic based on an analysis of variance, is computed. The pair with the highest probability Q(Flv a ,v b ) is merged. The 
parameters and limits of the merged line-cluster are recomputed using weighted linear regression. The whole cycle is 
repeated until no line-clusters pass the linear coherency test. 

30 

2.1.1 Linear Coherency 
Definitions: 

3S [0091] Linearly coherent: Two lines or line-clusters are said to be linearly coherent if they could potentially form part 
of the same side of the rectangular signal/shadow boundary. 

Linear agglomerative clustering: an iterative process for grouping together lines which are deemed (by some criteria) 
to be sufficiently similar. E.g. a pair-wise agglomerative clustering process based on an analysis of variance test, and 
broadly as described in (Silverman & Cooper), but with the mergable entities being lines rather than regions, and the 
40 constraint that two mergable regions be touching replaced by a constraint that two lines (or line-clusters) be linearly 
coherent. 

[0092] In the described embodiment of the present invention the linear coherency test is applied to each line or line- 
cluster pair during linear agglomerative clustering. For two line or line-clusters to be viable candidates for clustering: 

45 1 . They must be approximately collinear. 

2. Their extrapolations must overlap or touch. 

2.1.1.1 Collinearity 

so [0093] Two line-clusters are determined as collinear if: 

1 . The line-clusters are approximately parallel: 

55 Min (Q (1) WLR - Q (2) WLR ~ 2e7l) < 15 ° 

e in (-1.0.1) 
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and 

2. The dlstanc of closest approach betwe rt the two line-clusters Is less than some threshold: 



20 



25 



30 



35 



55 



d < 4 pixels. 



2.1.1.2 Overlapping lines 



[0094] Two line-clusters are said to have overlapping extrapolations if either has a perpendicular projection onto the 
10 other. 

2.1 .2 Analysis Of Variance 

[0095] The F-statistic based upon the analysis of variance is applied to test the likelihood that two sample distributions 
15 come from the same parent distribution. With reference to the definition of analysis of variance, and weighted linear 
regression, the chi -squared variates are the sum of square errors of the best fit straight line: 



where the line-cluster index j is any of a.b.m. The corresponding degrees of freedom are: 

v a = N a -2 



v b = N b -2 



where N a and N b are the numbers of edge-pixels in the lines (of labels a and b respectively). 
Steps iii) and iv) Construction and costing of hypotheses 



[0096] The role of the search tree is to find ail hypotheses with upto two line-clusters which could possibly be evidence 
for a multiply-exposed image. The reason for allowing upto two line-clusters for a multiple-exposure boundary is that 
the placement of X-ray masking material to protect each sub-image may not be perfectly aligned. When the masks are 
perfectly aligned, the result is a single edge as implied in figure 2. Figure 10 illustrates possible situations when the 
40 sub-image masks are not perfectly aligned. These correspond to hypotheses comprising two line-clusters and therefore 
at the third level of the search tree. 

[0097] A search tree is constructed of all legal combinations of line-clusters. The tree starts with a root node (no line- 
cluster) which corresponds to the null hypothesis that the image is a single-exposure, rather than multiple-exposure. 
[0098] The search tree has at most three levels, corresponding to the null hypothesis, and one level for each of the 
45 two line-clusters which can make up a multiple-exposure boundary. Each node of the search tree represents a hypoth- 
esis in its own right. The connection between the line-clusters, the hypotheses, and the nodes in the search tree are 
illustrated in figure 11. 

[0099] With reference to figure 1 1 , only the lines (a, b, c and d) are labelled which take part in hypothesis construction . 
In practice, an exhaustive search is made of all line-clusters, and line-cluster pair combinations, to test their possible 
so inclusion in a hypothesis. However, the search relies on a set of rules regarding placement of the line-clusters to select 
only those which could possible form part of a hypothesis. This includes, for example, rules to select only those clusters 
which are approximately central to the image, running approximately horizontally or vertically. 
[0100] Hypotheses are costed at the stage of construction. If there are n lines-clusters in the hypothesis, and if, for 
each cluster (index k), the fraction of its length for which there is edge data is f^ then the cost is : 



cost - -n-Max(l-f k ) wherein 0 < k < n. 0 <. n < 2 , 0 < f k jC l 
over k 
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[0101] For example : a null hypothesis has a default cost of 0; a hypothesis with on line-cluster with fraction 0.3 of 
its length supported by edg data has cost -1.7. Clearly, hypotheses are favoured which contain mor line-clusters, 
and for which there is more supporting edge data. 



s Step v) Evaluation of hypotheses 



[0102] Each hypothesis is propagated through the reasoning network. Each node of the network contains a rule 
which must return a true value in order that the hypothesis propagate onto the descendants of that node. Thus a rule 
value at a descendent node forms a logical AND with its antecedent. Rule values at nodes sharing a common ante- 
to cedent form a logical OR. An example is given in figure 8. Thus each hypothesis is subject to a complex logical ex- 
pression which must be true in order that the hypothesis is verified. 

[0103] The rule includes tests on the line-clusters and resulting sub-images comprising the hypothesis under eval- 
uation. These tests place constraints (upper and lower limits) on the following quantities : 



15 Geometry : Parallelism of line-cluster pairs (where there is more than one) 

Angle of line line-cluster(s) 

Centrality of line-clusters in the image 

Aspect ratio of resulting sub-images 
Line-cluster statistics Straightness 
20 Completeness (from one side of the image to the other) 

Contrast (normal to the cluster) 
Sub-image statistics Mean grey-level of mask overlap region (where applicable - see figure 10b) 

Mean grey-level of incomplete masking region (where 

applicable - see figure 10c). 



25 



35 



Step vi) Result selection 



[0104] All hypotheses which reach the terminal nodes of the reasoning network are deemed verified. Once all hy- 
potheses generated at step (iii) have been tested, that hypothesis with the least cost which has been verified at step 
30 (iv) is selected as the result. The default hypothesis is that the image is a single exposure only. This hypothesis is 
always verified as true, but has the highest possible cost (=0), and therefore is neglected in favour of any other feasible 
hypothesis. 



Step vii) Sub-image partitioning 

[0105] If the image is multiply-exposed, then the image is split into two sub-images consistent with the locus of the 
line-cluster(s) comprising the accepted hypothesis : 

If there is just one line-cluster, then this is used as a guide for the cut : if the line-cluster is exactly parallel to an 
40 image edge, then the cut follows its trajectory. If it is not exactly parallel to an image edge, then the line-cluster is 

bisected so as to ensure two perfectly rectangular sub-images. If there are two line-clusters, then a rectangular 
cut is made which otherwise is as central as possible in the masking overlap or incomplete masking bands. 
The result of this step is to create one (null hypothesis) or two rectangular sub-images which are individually passed 
on to the collimation detection process (steps vii to xi). In the case of the null hypothesis, it is to be understood (in 
45 the following) that the (single) sub-image is identically equal to the whole image. 

Steps viii) ix and x) Line detection and clustering 

[0106] Steps ix) and x) have already been described in detail (as steps i and ii) in EP-A-610 605. The steps are 
50 identical to the steps (i), and (ii) described above. They form the tow-level image processing and analysis applied to 
the sub-image. 

Steps xi) and xii) Construction and costing of hypotheses 

55 [0107] The role of the search tree is to find all hypotheses with uptofour line-clusters which could possibly be evidence 
for a signal/shadow boundary. 

[0108] A search tree is constructed of all legal combinations of line-clusters. The tree starts with a root node (no line- 
cluster) which corresponds to the null hypothesis that there is no signal/shadow boundary, and therefore no collimation, 
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in the current sub-image. 

[0109] The search tre has at most five levels, corresponding to the null hypothesis, and one level for each of th 
four line-clusters which can make up a signal/shadow boundary. Each node of the search tr e represents a hypoth sis 
in its own right. The connection between the line-clusters, the hypotheses, and the nodes in the search tree are illus- 
5 trated in figure 12. 

[0110] With reference to figure 12 an exhaustive search is made of all line-clusters, and line-cluster combinations 
(upto four deep), to test their possible inclusion in a hypothesis. However, the search relies on a set of rules regarding 
the geometry of the line-cluster combination to select only those combinations with could possibly form part of a hy- 
pothesis. This includes, for example 

10 

a) Rules to select only those clusters which are mutually orthogonal (sides of a rectangle). 

b) Rules to select combinations which form a closed boundary either between themselves, or including the edge 
of the sub-image. 



1S [0111] It is through the enforcement of the closure rules that, for instance, none of clusters b, c, d, g, or e may form 
a one-sided hypothesis. 

[0112] Hypotheses are costed at the stage of construction. If there are n line-clusters in the hypothesis, and if, for 
each cluster (index k), the fraction of its length for which there is edge-data is f^ then the cost is : 

20 

cost - -n-Max(l-fk) wherein 0 <, k < n, 0 < n < 4, 0 < < 1 
over k 



[0113] For example : a null hypothesis has a default cost of 0; a hypothesis with one line-cluster with fraction 0.3 of 
25 its length supported by edge data has cost -1.7. Clearly, hypotheses are favoured which contain more line-clusters, 
and for which there is more supporting edge data. 



Step xiii) Evaluation of hypotheses 



30 [0114] Each hypothesis is propagated through the reasoning network. Each node of the network contains a rule 
which must return a true value in order that the hypothesis propagate onto the descendants of that node. Thus a rule 
value at a descendent node forms a logical AND with its antecedent. Rule values at nodes sharing a common ante- 
cedent form a logical OR. An example is given in figure 8. Thus each hypothesis is subject to a complex logical ex- 
pression which must be true in order that the hypothesis is verified. 

35 [0115] The rule includes tests on the line-clusters and resulting sub-images comprising the hypothesis under eval- 
uation. These tests place constraints (upper and lower limits) on the following quantities : 



40 



so 



Geometry : 



Line-cluster statistics 



45 Region statistics 



Step xiv) Result selection 



Orthogonality of line-clusters (when there is more than one) Orientation of the line- 
clusters) 

CentraYity of resulting signal region in the sub-image 
Aspect ratio of resulting signal region 
Straightness 

Completeness (from one side of the image to the other) 

Contrast (normal to the cluster) 

Mean grey-level of signal and shadow regions 

Variance of signal and shadow regions 

Edge-point density of signal and shadow regions 

Mean grey-level of half planes on shadow side of each cluster 

Edge-point density in half planes on shadow side of each cluster. 



[0116] All hypotheses which reach the terminal nodes of the reasoning network are deemed verified. Once all hy- 
potheses generated at step (xi) have been tested, that hypothesis with the least cost which has been verified at step 
55 (xii) is selected as the result. The default hypothesis is that there is no collimation, and that th current sub-image is 
therefore completely exposed. This hypothesis is always verified as true, but has the high st possible cost (=0), and 
therefore is neglected in favour of any other feasible hypothesis. 
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Step xv) Binary mask image 

[0117] Th binary mask is an image constructed by setting to some int ger value all pix ts in a signal region, as 
determined by successive iterations of steps vii) to xi). The sam integer is assigned to all pixels In a signal region; a 
s different integer for each signal region. The remaining (zero-valued) pixels are deemed to be either in the shadow 
region(s), in a multiple-exposure mask overlap region (figure 10b), or in a multiple exposure region of incomplete 
masking (figure 10c). This mask image is then presented to the histogram equalisation and image display stages to 
act as a mask for further processing and display functions (to be executed upon the original X-ray image). 

10 Network Property a) Archetype Classification 

[0118] The multiple-exposure and signal/shadow reasoning networks each embody an archetype classification and 
refinement philosophy. By archetype classification we mean that the higher (earlier) branches of the network are de- 
signed to channel (propagate) the current hypothesis according to its classification as one of a set of well defined 

15 topological and geometrical configurations. In the case of the multiple -exposure network, these archetypes are illus- 
trated in figure 1 3. For the case of the signal/shadow reasoning network, these archetypes are illustrated in figure 1 4. 
[0119] The rules at the nodes in the higher branches still belong to the sets described in steps iv and x respectively. 
However, the role of these rules is to act as a filter with the aim of having the current hypothesis proceed to the deeper, 
lower nodes in just one of the major branches. In that sense, the current hypothesis is deemed to be classified according 

20 to an archetype. 

Network Property b) Refinement 

[0120] Once a hypothesis has been classified, further rules must be satisfied before it is deemed verified. These 
25 refinement rules seek to exclude all those hypotheses which, for the given topological and geometrical classification, 
broadly speaking, do not have the requisite line-cluster or region statistics. These statistical constraints are those which 
have been found empirically, from a network training phase, to be associated with true hypotheses. Further geometrical 
rules may also be applied during refinement where appropriate. 

30 Network Property c) Training 

[0121] The thresholds for the rules at all stages in the network are obtained from interactive supervision of the pas- 
sage of good and bad (true and false) hypotheses through the network. Changes are made to the rule thresholds 
consistent with the intended function of rules to classify or refine/reject hypotheses. The aim of this training is to arrive 
3S at a network for which all true hypotheses, and more of the false hypotheses, survive the refinement stages. The 
network contains a memory of previously supervised hypotheses, so that adjustments during supervision are consistent 
with previous training passes, and the network converges towards a stable behaviour. 



40 Claims 

1. A method of determining a multi-exposure boundary between sub-images in a multi-exposure X-ray image repre- 
sented by a digital signal representation and of determining a signal/shadow boundary in each sub-image, com- 
prising the steps of 

45 

i) Extracting low-level line primitives from the X-ray image, 

ii) Forming intermediate-level primitives from the low-level primitives, 

iii) Building hypotheses as to the location of the multiple-exposure boundary from combinations of intermediate- 
level primitives, 

so jv) Associating with each hypothesis a cost, 

v) Subjecting each hypothesis to evaluation tests via a rule-based partitioning reasoning network, 

vi) Selecting the hypothesis with the least cost, 

vii) Partitioning the image into sub-images according to the selected hypothesis, and 

viii) Submitting each sub-image for subsequent detection of the signal/shadow boundary to the following steps 
55 (be) to (xv): 

ix) Extracting low-level primitives from each sub-image, 

x) Forming intermediate-level primitives from the extracted tow-level primitives, 

xi) Building hypotheses as to the location of the signal/shadow boundary from combinations of intermediate- 
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level primitives, 

xii) associating a cost with each hypothesis, 

xiii) Subjecting each hypoth sis to evaluation tests via a rule bas d collimation reasoning network, 

xiv) . Selecting a hypothesis with th least cost, 

* xv) Construction of a binary-valued mask image representing signal and shadow area to be used to extract 

the signal area in the X-ray image. 

A method according to claim 1 , wherein said low level primitives are lines, said intermediate level primitives are 
line-groups and wherein said hypotheses for the multiple-exposure boundary are formed from every possible line- 
cluster and pair of line-clusters, including, in addition, the hypothesis that there is no multiple-exposure boundary. 

A method according to claim 1 , wherein the evaluation tests applied to a hypothesis as to the location of the 
multiple-exposure boundary, comprise any logical combinations of the following tests for : 

- boundary closure, for which test to be passed, line-clusters in the hypothesis under consideration, when having 
been extrapolated must each terminate on the edge of the image, 

boundary location, for which test to be passed, said line-clusters must pass close to the centre of the image, 
boundary alignment, for which test to be passed, said line-clusters must be nearly parallel to a side of the image, 
boundary completion, for which test to be passed, the fraction of the boundary perimeter length, for which 
there exists evidence from the low-level primitives, must be greater than some threshold. 

A method according to claim 2 wherein said cost that assigned to each multiple-exposure hypothesis is lowest if 
a hypothesis is that of two approximately parallel line-clusters, and wherein said cost increases as the number of 
visible line-clusters in a hypothesis decreases, and the cost is highest if the hypothesis has no line-clusters. 

A method according to claim 1 , wherein the first stages of evaluation involve detection of topological and geomet- 
rical archetypes for the multiple-exposure, and signal/shadow boundaries. 

A method according to any of the preceding claims wherein image processing performed on said digital signal 
representation of said X-ray image is limited to pixel values of said image within the signal/shadow boundary. 

A method according to claim 1 wherein a digital signal representation of an X-ray image is obtained by scanning 
an exposed photostimulable phosphor screen with stimulating irradiation, detecting light emitted upon stimulation 
and converting the detected light into a digital signal representation. 



Patentanspruche 

1 . Ein Verfahren zum Bestimmen einer Mehrf achbelichtungsgrenze zwischen Teilbildern in einem durch eine digitale 
40 Signaldarstellung dargestellten Mehrfachbelichtungs-Rontgenbild und zum Bestimmen einer Signal-Schatten- 

Grehze in jedem Teilbild, wobei das Verfahren folgende Schritte umfaBt: 

i) Ableiten von Liniengrundelementen niedriger Stufe aus dem Rdntgenbild, 

ii) Bilden von Grundelementen der Zwischenstufe aus den Grundelementen niedriger Stufe, 

45 jjj) Bilden von Hypothesen hinsichtfich des Ortes der Mehrfachbelichtungsgrenze aus Kombinationen von 

Grundelementen der Zwischenstufe, 

iv) jeder Hypothese Kosten zuordnen, 

v) jede Hypothese uber ein Unterteilungsargumentationsnetz auf Regelbasis Bewertungspruf ungen unterzie- 
hen, 

so vi) Auswahl der Hypothese mit den geringsten Kosten, 

vii) Unterteilen des Bilds in Teilbilder gemaB der ausgewahlten Hypothese, und 

viii) Unterziehen jedes Teilbikte zur spateren Erfassung der Signal-Schatten-Grenze den foJgenden Schritten 
(ix) bis (xv): 

ix) Ablerten von Grundelementen der niedrigen Stufe aus jedem Teilbild, 

ss x ) Bilden von Grundelementen der Zwischenstufe aus den abgeleiteten Grundelementen niedriger Stufe, 

xi) Bilden von Hypothesen hinsichtlich des Ortes der Signal-Schatten-Grenze aus Kombinationen von Grund- 
elementen der Zwischenstufe, 

xii) jeder Hypothese Kosten zuordnen, 
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xiii) jed Hypotheso Ober ein Kollimationsargumentationsnetz auf Regelbasis Bew rtungsprOfungen unterzie- 
hen, 

xiv) Auswahl einer Hypothes mit den geringsten Kosten, 

xv) Konstruier n eines binarw rtigen Maskenbilds, das Signal- und Schattenbereich darstelft, zur vferwendung 
5 zum Ablerten des Signaibereichs in dem Rontgenbitd. 

2. Verfahren nach Anspruch 1 , wobei es sich bei den Grundelementen niedriger Stufe um Linien handelt, die Grund- 
elemente der Zwischenstuf e Liniengrupen sind und wobei die Hypothesen f Or die Mehrfachbelichtungsgrenze aus 
jedem moglfchen LinienbOndel und Paar von LfnienbOndeln gebildet werden, einschlieBlich auBerdem der Hypo- 

10 these, daB keine Mehrfachbelichtungsgrenze vorliegt. 

3. Verfahren nach Anpruch 1 , bei dem die auf eine Hypothese angewendeten Bewertungsprufungen hinsichtlich des 
Orts der Mehrfachbelichtungsgrenze alle beliebigen logischen Kombinationen der folgenden PrOf ungen auf fol- 
gende Punkte umfassen: 

75 

GrenzschlieBung, wobei um diese PrOfung zu bestehen, die LinienbOndel in der betrachteten Hypothese, 
wenn sie extrapoliert worden sind, jeweils am Rand des Bildes enden mOssen, 
- Grenzort, wobei um diese Prufung zu bestehen, die LinienbOndel nahe an der Mitte des Bilds vorbeigehen 
mussen, 

20 - Grenzausrichtung, wobei um diese PrOfung zu bestehen, die LinienbOndel fast parallel zu einer Seite des 
Bildes veriaufen mussen, 

Grenzvollendung, wobei um diese Prufung zu bestehen, der Anteil der Grenzumfangslange, fOr den Befund 
von den Grundelementen niedriger Stufe vorliegt, Ober einem gewissen Schwellwert liegen muB. 

25 4. Verfahren nach Anspruch 2, wobei die jeder Mehrfachbelichtungshypothese zugewiesenen Kosten am niedrigsten 
sind, wenn es sich bei einer Hypothese um eine von zwei ungefahr parallelen LinienbOndeln handelt, und wobei 
die Kosten mit der Abnahme der Anzahl sichtbarer LinienbOndel in einer Hypothese ansteigen und die Kosten am 
hochsten sind, wenn die Hypothese keine LinienbOndel aufweist. 

30 5. Verfahren nach Anspruch t , bei dem bei den ersten Stadien der Bewertung die topologischen und geometrischen 
Archetypen fOr Mehrfachbelichtung und Signal-Schatten-Grenzen erfaBt werden. 

6. Verfahren nach einem der vorhergehenden AnsprOche, bei dem die Bildverarbeitung, die an' der digitalen Signal- 
darstellung des Rontgenbildes durchgefuhrt wird, auf Bildpunktwerte des Bildes innemalb der Signal-Schatten- 

35 Grenze begrenzt ist. 

7. Verfahren nach Anspruch 1, bei dem eine digitate Signaldarstellung eines Rontgenbildes durch Abtasten eines 
belichteten photostimulierbaren Leuchtstoffschirms mit stimulierender Strahlung, Erfassen von bei Stimulierung 
emittiertem Licht und Umwandeln des erfaBten Lichts in eine digitale Signaldarstellung erhalten wird. 

40 

Revendications 

1. Proc6d6 de determination cfune fronti^re cfexposition multiple entre des sous-images dans une image radiogra- 
45 phique k exposition multiple representee par une representation de signal num6rique, et de determination d'une 

frontidre signal/ombre dans chaque sous-image, comprenant les etapes consistant k: 

i) Extraire des primitives de lignes de bas niveau de I'image radiographique, 

ii) Former des primitives de niveau intermediate & partir des primitives de bas niveau, 

so jjj) Elaborer des hypotheses quant & Pemplacement de la frontiere d'exposition multiple h partir de combinai- 

sons de primitives de niveau intermediate, 

iv) Associer chaque hypothese a un cout, 

v) Soumettre chaque hypothese a des tests devaluation par le biais d'un reseau de raisonnement par partage 
k base de regies, 

55 vi) Seiectionner Thypothese de moindre cout, 

vii) Partager I'image en sous-images selon Phypothese seiectionnee, et 

viii) Soumettre chaque sous-image aux etapes (ix) k (xv) suivantes en vue de la detection subs6quente de la 
fronti6re signal/ombre : 
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ix) Extraire des primitives de bas niveau de chaque sous-image, 

x) Former des primitives de niv au intermediate k parti r des primitives de bas niveau xtrait s, 

xi) Elaborer des hypotheses quant k Pemplacement de la frontier signal/ombr 6 partir des combinaisons de 
primitives d niveau intermediate, 

xii) Associer un cout k chaque hypothese, 

xiii) Soumettre chaque hypothese k des tests devaluation par !e biais d'un r6seau de raisonnement de colli- 
mation k base de regies, 

xiv) S6lectionner une hypothese de moindre coOt, 

xv) Construire une image de masque k valeurs binaires representant ta zone de signal et d'ombre k utiliser 
pour extraire ta zone de signal dans I'image radiographique. 

Precede selon la revendication 1 , dans lequel lesdites primitives de bas niveau sont des lignes, lesdites primitives 
de niveau intermediate sont des groupes de lignes, et dans lequel lesdites hypotheses pour la frontiere d'exposition 
multiple sont form6es k partir de chaque regroupement de lignes et paire de regroupements de lignes possibles, 
y compris, en outre, I'hypothese qu'il n*y a pas de frontiere ^exposition multiple. 

Proc6d6 selon la revendication 1, dans lequel ies tests devaluation appliques k une hypothese quant k I'empla- 
cement de la fronti6re d'exposrtion multiple comprennent toute combinaison fogique des tests suivants de: 

fermeture de la frontiere; pour qu'un tel test sort passe avec succes, Ies regroupements de lignes dans I'hy- 
pothese consideree doivent, une fois extrapoies, chacun aboutir sur le bord de I'image, 
emplacement de la frontiere; pour qu'un tel test sort passe avec succes, lesdits regroupements de lignes 
doivent passer k proximite du centre de I'image, 

alignement de la frontiere; pour qu'un tel test soit passe avec succes, lesdits regroupements de lignes doivent 
etre pratiquement paralieies k un cdte de i'image, 

achevement de la frontiere; pour qu'un tel test soit passe avec succes, la fraction de la longueur du p6rimetre 
de la frontiere, pour laquelle il existe des preuves provenant des primitives de bas niveau, doit etre sup6rieure 
k un certain seuil. 

Proc6d6 selon la revendication 2, dans lequel ledit cout affecte k chaque hypothese d'exposition multiple est le 
plus bas si une hypothese est cede de deux regroupements de lignes approximatrvement paralieies, et dans lequel 
ledit cout augmente k mesure que le nombre de regroupements de lignes visibles dans une hypothese diminue, 
et le cout est le plus eieve si I'hypothese n'a aucun regroupement de lignes. 

Procede selon la revendication 1, dans lequel Ies premiers stades de revaluation mettent en jeu la detection 
d'archetypes topologiques et geom6triques pour I'exposition multiple et Ies frontiers signal/ombre. 

Proc6d6 selon I'une quelconque des revendications pr6c6dentes, dans lequel le traitement de I'image realise sur 
ladrte representation de signal numerique de iadite image radiographique est limite aux valeurs de pixels de ladite 
image k I'int6rieur de la frontiere signal/ombre. 

Procede selon la revendication 1 , dans lequel une representation de signal numerique d'une image radiographique 
est obtenue par balayaged'un ecran luminescent photostimulable expose a Paide d'un rayonnement de stimulation, 
detection de la lumiere emise suite k la stimulation, et conversion de la lumiere detectee en une representation 
de signal numerique. 
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